How is the magnetic reconnection derived from magnetohydrodynamics equations? 
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We clarify how magnetic reconnection can be derived from magnetohydrodynamics (MHD) equations in a 
way that is easily understandable to university students. The essential mechanism governing the time evolution 
of the magnetic field is diffusion dynamics. The magnetic field is represented by two components. It is clarified 
that the diffusion of a component causes a generation of another component that is initially zero and, accordingly, 
that the magnetic force lines are reconnected. For this reconnection to occur correctly, the initial magnetic field 
must be directed oppositely in the two regions, e.g., y > and y < 0; must be concave (convex) for y > 
(y < 0); and must be saturated for y far from the x axis, which would indicate the existence of the current sheet. 
It will be clear that our comprehension based on diffusion runs parallel to the common qualitative explanation 
about the magnetic reconnection. 



I. INTRODUCTION 

A solar flare is the most violent explosion in the solar sys- 
tem. It occurs over the Sun's surface, resulting in the bright- 
ening of electromagnetic waves over a wide wavelength range 
and the release of about 10 22 -10 25 joule of energy within 
about 10 2 -10 3 s|l. 

How such a large amount of energy is released in such a 
short time has puzzled scientists for many years. One of the 
epoch-making approaches to this problem is to exchange the 
magnetic energy into kinetic energy and the heat of plasma by 
reconnecting magnetic field lines. This is referred to as mag- 
netic reconnection. The time evolution of the magnetic field 
is determined by a magnetic induction equation that includes 
two contrary effects, "freezing" and diffusion of the magnetic 
field, which is explained in detail later. Diffusion of the mag- 
netic field is caused by Joule dissipation. If the electric resis- 
tivity is low, the magnetic field is "frozen" into the plasma. 
However, if an anomalous resistivity is induced for some rea- 
son and the diffusion effect becomes dominant, the freezing of 
the magnetic field ceases and its structure can be altered. 

In other words, the diffusion causes structural change of 
the magnetic field. However, the time for the release of 
the magnetic energy is quite long. In the case of the solar 
flare, the diffusion constant rj is about 1 m 2 /s and the typi- 
cal scale of length for the energy release, which we shall de- 
note by Lo, is 10 8 m. Therefore, the time can be estimated 
as Lq 2 It] ~ 10 16 s ~ 10 8 yr, which is not comparable with 
the observed release time[2]. Hence, we can explain the heat- 
ing of plasma and the release of a large amount of energy by 
considering only the time evolution of the magnetic field, but 
cannot explain the release time. 

For resolving this problem, Dungey proposed a model con- 
sidering the electromagnetic field caused by the motion of 
the plasma and, thus, using magnetohydrodynamics (MHD) 
with negligible gas pressure. He showed that a current sheet 
forms between opposite directions of the magnetic fieldQ |5fl. 
Dungey was also the first to suggest that the configuration 
of the magnetic field can change. Subsequently, Sweet and 



Parker expanded on Dungey's model by considering the in/out 
flow of the plasma and using MHD equations; this model is re- 
ferred to as Sweet-Parker's model^^^^. According to their 
model, a plasma flowing into the current sheet with a velocity 
v, flows out along the sheet with the Alfven velocity Va ■ The 
velocity v, satisfies the following relation. 
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where R m denotes the magnetic Reynolds number, which can 
be expressed as 



Rm — 



LoV A 



(2) 



The magnetic energy moves with velocity v,, and the time with 
which the energy is released over a length Lq is given by 



(3) 
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If we use 10 8 m and 1 m 2 /s as values for Lq and rj, respec- 
tively, and estimate the Alfven velocity as 10 8 m/s, the energy 
release time becomes 10 8 s. Therefore, we can find that the 
time is much shorter than that for the model considering only 
diffusion dynamics. However, this time is still longer than 
what is observed. Petschek reduced the reconnection region 
in order to resolve this discrepancy, resulting in the release 
time of 10 2 -10 3 s, which is closer to the observed oneH. 

As mentioned above, the theory of magnetic reconnection 
has been developed. However, the concept used in describ- 
ing the time evolution of the magnetic field is that of diffu- 
sion. Hence, the question arises as to why diffusion produces a 
structural change of the magnetic field. Of course, we can of- 
ten find a qualitative explanation by drawing analogy between 
magnetic field lines and rubberbands: the energy of short con- 
nected magnetic field lines is lower than that of extended mag- 
netic field lines, as is shown in Fig.[TJ However, this analogy 
does not account for the aspect of diffusion dynamics. One 
of ways to determine that magnetic reconnection is caused by 
diffusion is to calculate the magnetic induction equation and 
other MHD equations numerically. However, such a calcula- 
tion would be required to be done over a complex computer 
program using an expensive, high-specification machine. In 
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FIG. 1: Illustration of magnetic reconnection by using an analogy 
between magnetic field lines and rubber bands. Enormous energy is 
stored in the extended magnetic field lines shown in (a). The occur- 
rence of dissipation makes the magnetic force lines reconnect, as in 
(b), resulting in the release of stored energy. 



this paper, we demonstrate how the magnetic field reconnects 
not through such an advanced technical process, but through 
a method that can be easily understood by university students 
who are familiar with diffusion phenomena. 

This paper is organized as follows. First, we derive the 
equations governing the time evolution of magnetic and ve- 
locity fields and define their initial and boundary conditions 
in Sec. [II] Next, the diffusion phenomena are revisited in 
Sec. IIH Al after which we demonstrate how the magnetic field 
lines are reconnected using the basic knowledge of diffusion 
in Sec. IIIIBl Finally, in Sec. lIVl we conclude this paper and 
indicate the equality of our comprehension and the common 
qualitative explanation. 



II. FUNDAMENTAL EQUATIONS 

In this section, we enumerate the equations governing the 
time evolution of magnetic and velocity fields. For simplicity, 
these fields are assumed to exist only on the x-y plane. 

The time evolution of the magnetic field, B, is expressed by 
the following magnetic induction equation. 



dB_ 



Vx(vxB) + t]AB 



(4) 



where v is the velocity field of the plasma and rj is the mag- 
netic diffusivity, which can be defined by using the electric 
resistivity, rj e , and the magnetic permeability, [i, as r\ = r\ e ln. 
The right hand side has two terms: the first represents the 
effect of the freezing of the magnetic field into the plasma, 
and the second represents the diffusion. For the case with 
r\ — 0, i.e., the case of a perfect conductor, the magnetic field 
moves together with the plasmaQlIS]. This magnetic induc- 



tion equation can be derived by incorporating Faraday's law 

-V x E , (5) 
the electric field of the Ohm's law 



dB 



and Ampere's law 



T] e j = E + vxB 



VxB=nj, 



(6) 



(7) 



in which j is the current density. See Appendix for the details. 

Next, the equation of motion of the plasma that is consid- 
ered as a perfect incompressible fluid is 

Pm|J- t -(v-V)vj = -V/> + /XB, (8) 

where p m and p represent the mass density and the pressure 
of plasma, respectively. The above two equations, Eqs and 
©, are sufficient for us to understand the magnetic reconnec- 
tion. 

We show the change in the electric resistivity, r\ e , with 
time, which is how the structural change of the magnetic 
field starts, as follows: before t — 0, the value is zero ev- 
erywhere; subsequently, it increases drastically only over the 
region —L < x < L and —w < y < w, which we call the recon- 
nection region. Hence, the magnetic field "frozen" into the 
plasma starts to "melt" at t — 0. Then, we set the change in 
magnetic diffusivity, rj, in this region as follows: 



n = 



(f < 0) 
//* (» 1) (r > 0) 



(9) 



In addition, we shall call the other region, \x\ > L or |y| > w, 
the frozen region. 

The initial magnetic field, B(r, t < 0), is chosen to be par- 
allel with the x axis at y > and antiparallel with it at y < 0. 
Thus, we adopt an odd function as the magnetic field. Here, 
we use the following simple expression 

B(r, t < 0) = (B x (r, t < 0), B y (r, t < 0), 0) 

= B (tanhay,0,0), (10) 

where a is some constant whose dimension is the inverse of 
length. 

From Ampere's law, the current density can be obtained as 

j(r,t < 0) = -VxB(r,t< 0) 

= -(0, 0, 8 x B y (r, t < 0) - d y B x (r, t < 0)) 
H 

= _^(0,0,sech 2 ay) . (11) 

This current density localizes around the x axis, that is, the 
current sheet. The existence of this current sheet between op- 
positely directed magnetic fields is important in the magnetic 



reconnection. Therefore, we choose the hyperbolic tangent 
among odd functions: the derivative of the initial B x with re- 
spect to y must localize around the x axis. In other words, 
this initial magnetic field must be concave (convex) at y > 
(y < 0) and be saturated for y far from the x axis. 

In this paper, we regard the plasma as a perfect conductor 
with no resistivity until t — 0, i.e., rj — r\ e — 0. The velocity 
field must be zero in order to keep this magnetic field static: 

v(r, t < 0) = (v x (r, t < 0), v y (r, t < 0), v z {r, t < 0)) 

= (0,0,0). (12) 

This is because the plasma and the magnetic field move to- 
gether for the case with 77 = 0. 

For the velocity field to maintain this state, the following 
relation, obtained from Eq. (|8), must be satisfied. 

-Vp = -jxB (13) 

From Eq. (fTTT i. the right hand side of Eq. (Q~3) becomes 

-j(r, t < 0) x B(r, t < 0) 

= (j z (r, t < 0)B,(r, t < 0), -j z (r, t < 0)B x (r, t < 0), 0) 
= — — (0, sech 2 ay tanhay, 0) . (14) 

Hence, we obtain the pressure gradient, 
aB 2 

- Vp = — - (0, sech 2 ay tanh ay, 0) . (15) 

These magnetic and velocity fields do not change in the frozen 
region after t — 0. 



III. TIME EVOLUTION OF MAGNETIC AND VELOCITY 
FIELDS 

A. diffusion 

In this subsection, we revisit the diffusion dynamics. Al- 
though diffusion is not the primary topic of this paper, the 
main equation governing the time evolution of the magnetic 
field is the diffusion equation, as shown in Eq. 

The diffusion equation in one dimension is 

d l P(x,t) = Dd 2 t P(x,t), (16) 

in which D is the diffusion constant. We can interpret this 
equation as follows: P(x, f) increases (decreases) with time 
if the second derivative with respect to the spatial variable 
is positive (negative). P(x, t) is convex (concave) at a region 
where the second derivative is positive (negative). Therefore, 
a picture of the diffusion dynamics is a rise of the convex re- 
gion and a fall of the concave one, as is shown in Fig. [2] 
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FIG. 2: Time evolution of the diffusion equation, Eq. l !16t . from an 
initial state (gray line) to a later state (black line). The concave region 
of P(x, t) decreases, and the convex one swells up. 



B. magnetic reconnection 

Here, we investigate how the magnetic field changes for a 
small t by the occurrence of dissipation at t — 0. 

From the values of the magnetic and velocity fields at t < 0, 
shown as in the previous section, one can obtain 

d t B\ t=0 = tfAB (17) 

and 

d t v\ t=0 = (18) 

From these equations, it is clarified that the magnetic field 
evolves with time earlier than the velocity field. Hence, let us 
consider the time evolution of the magnetic field maintaining 
the value of the velocity field as that in Eq. (fT2l) . and then, let 
us investigate how this time evolution affects the zero velocity 
field. 

Because we can regard B v at a small t as being nearly equal 
to zero, magnetic induction equations at the reconnection re- 
gion become 

d,B x (r, t) = -dyv y (r, t)B x {r, t) + rfhB x {r, t) (19) 

and 

d t B y (r,t) = d x v y (r,t)B x (r,t). (20) 

Moreover, from Eqs. (TTOt and (Q~2), we can regard B x and v, 
at a small t as 

B x (r, t) =i Bq tanh ay (21) 

and 

v y (r, - . (22) 

Therefore, the right hand side of Eq. (l20t can be considered to 
be zero at such a time stage. That is to say, B x changes with 
time earliest among the components of the magnetic field. 

By solving Eq. (T% , we obtain the behavior of B x only at 
y > 0, as the initial B x is odd under the transformation y — > -y 
and this equation is invariant under transformations B x — > ~B X 
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and y — > —y. 111 111 Therefore, the relation B x (x, 0, t) = is 
always satisfied. 

Equations (TT~9b and (T22l i together compromise the exact dif- 
fusion equation. Because the initial magnetic field, B x (r, 0) = 
Bo tanh ay, is concave for y > 0, the effect that results in dent- 
ing of the convex region works in the same way as that ex- 
plained in Sec. IIII Al As mentioned in the above paragraph, 
however, B x is always zero at y = 0. Moreover, the value 
of the magnetic field is fixed at y = w and x - +L, as these 
are boundaries between the reconnection and frozen regions. 
Therefore, the denting is remarkable at the center of the re- 
connection region with y > 0, as shown in Fig.[3j this denting 
is indicated by 



d x B x 



>0 (forx>0) 
<0 (forx<0) 



(23) 



This absolute value increases with the distance from the y axis, 
which can be seen easily from Fig. [3J a )- 



(a) 



B x 
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^B x (x,y,t=0) 
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x B x (x,y,t>0) 
>- 




>■ y 



of the magnetic field as 

Ur,t) = -{d x B y (r,t)-d y B x ( r ,t)} , 



(24) 



changes with time in the following way: d y B x becomes 
smaller than its initial value, i.e., aBo sech 2 ay, particularly 
at a small y because of the denting effect of diffusion (See 
Fig. [3jb)). Therefore, the absolute value of j z localizing in 
the vicinity of the x axis at t — diminishes. The follow- 
ing qualitative explanation may help us understand this result 
easily: The occurrence of dissipation reduces the current den- 
sity. Hence, the balance between the pressure gradient and 
the force from the magnetic field described by Eq. (TPTt breaks 
down; the time evolution of the velocity field is thus governed 
only by the pressure gradient. Because this gradient is shown 
as Eq. ( TTBT ). only the y component of the velocity field has the 
following value, which does not depend on x: 



Vy 



>0 (fory>0) 
<0 (fory<0) 



(25) 



This absolute value becomes smaller with the distance from 
the x axis. The reason is as follows. The absolute value of the 
current density is originally small at a region far from the x 
axis, and hence, the pressure gradient balanced with the force 
caused by the current density and the magnetic field is also 
small at this region. 

As we have seen, d x B x and v y have non-zero values after 
t — owing to the diffusion dynamics in spite of the fact 
that these initial values are zero. From Eq. (T20t . the change 
of By with time is determined by the product of d x B x and v y , 
since v y is independent of x. In short, B Y starts to evolve with 
time since the right hand side of Eq. (l20b is no longer zero. 
Considering B v (r, 0) = 0, Eq. d23l , and v y > for y > 0, one 
can reach the following conclusion: 



B y 



>0 (forx>0) 
< (for x < 0) 



(26) 



This absolute value increases with distance from the y axis. 
On the other hand, this value decreases with distance from the 
x axis because of the behavior of v y as mentioned before. 

Hence, we can understand that the magnetic field directed 
parallelly or antiparallelly along the x axis at t — starts to 
reconnect owing to the occurrence of B y , which is described 
in Fig. g] 



FIG. 3: (color online) Time evolution of B x and its derivative with 
respect to y in the reconnection region. The red and black curves 
respectively represent B x (x, y, t > 0) and S,(x, y,t = 0), that is, 
Bo tanh ay. The dashed lines in (b) are the derivatives of the curves 
at y = 0. From both figures, it is clarified that a denting effect of the 
diffusion process works. In particular, we can find from (b) that the 
slope of the red dashed curve denoting d y B x (x,y, t = 0) is larger than 
that of the black curve denoting d y B x (x,y, t > 0). 

Because of the time evolution of B x , the z component of the 
current density, j z , which is described by using the derivatives 



IV. CONCLUDING REMARKS 

We have clarified how the diffusion gives rise to the mag- 
netic reconnection. The magnetic induction equation has two 
components. First, the x component, B x , diffuses, and then, 
the y component, B y , which is initially zero, is generated by 
the time evolution of B x . For the magnetic field lines to be 
reconnected correctly, the initial B x must be concave (convex) 
at y > (y < 0) and must be saturated for y far from the x axis 
to such an extent as 1 /a as in this paper; this behavior would 
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FIG. 4: Illustration of magnetic reconnection based on our com- 
prehension: (a) initial magnetic field described by Eq. d 1 0> where 
By = 0, (b) diffusion of B x causes emergence of B y , denoted by gray 
arrows, and (c) the composed magnetic field. 



pation and magnetic reconnection is generated. The common 
qualitative explanation of the magnetic reconnection using the 
current sheet is connected to our comprehension based on dif- 
fusion in this paper. Finally, we emphasize that if we adopt 
other concave (convex) functions of y at y > (y < 0), which 
are saturated as the initial magnetic field, magnetic reconnec- 
tion can occur. 



Appendix A: Magnetic induction equation 

By substituting the electric field of Ohm's law, 
rjej = E + v x B , 
into the left-hand side of Faraday's law, 



we can obtain 



— = -VxE, 

dt 



— = Vx(vxfi)-// t Vxj. 

at 



Ampere's law, i.e., 



V x B = fij 



(Al) 



(A2) 



(A3) 



(A4) 



relates the current density to the magnetic field. Therefore, 
equation dA3l ) becomes 

— = V x (v x B) - — V x (V x B) . (A5) 

Of fi 

For vector a, the following relation is satisfied. 

V x (V x a) = V(V ■ a) - Aa (A6) 

With regard to the magnetic field, V • B = is satisfied, as 
there is no magnetic monopole. Thus, 

V x (V x B) = -AB , (A7) 

and so, the induction equation is obtained as follows: 

dB n e 

— = V x (v x B) + — aB . (A8) 

at /i 



indicate that the current density is localized around the x axis, 
i.e., the current sheet. This current sheet is also necessary in 
the qualitative explanation. If an anomalous resistivity occurs, 
the current dissipates and the magnetic field starts to melt. As 
a result, enormous energy breaks out through the Joule dissi- 



Acknowledgments 

We would like to thank Hiroe Minagawa, Prof. Morikawa 
Masahiro, and members of Astrophysics Laboratory at 
Ochanomizu University for the extensive discussions. 



[1] T. Yokoyama, "Simulations of the magnetic reconnection as a [2] E. R. Priest, and T. G. Forbes, Magnetic Reconnection: MHD 
basic process of the solar activities", Nagare 29, 33-38 (2010). Theory and Applications (Cambridge University Press, Cam- 



6 



bridge, 2000). 

[3] R. G. Giovanelli, "A theory of chromospheric flares", Nature 
158, 81-82 (1946). 

[4] F. Hoyle, Some Recent Researches in Solar Physics (Cam- 
bridge University Press, Cambridge, 1949). 

[5] J. W. Dungey, "Conditions for the occurrence of electrical dis- 
charges in astrophysical systems", Phil. Mag. 44, 725-738 
(1953). 

[6] P. A. Sweet, "The neutral point theory of solar flares", Electro- 
magnetic Phenomena in Cosmical Physics, IAU Symp. 6, 123- 
134 (Cambridge University Press, Cambridge, 1958). 



[7] P. A. Sweet, "The production of high-energy particles in solar 

flares", Nuovo Cimento Suppl. 8 Ser. X, 188-196 (1958) . 
[8] E. N. Parker, "Sweet's mechanism for merging magnetic fields 

in conducting fluids", J. Geophys. Rev. 62, 509-520 (1957). 
[9] H. E. Petscheck, "Magnetic field annihilation", Physics of Solar 

Flares, ed. W. N. Hess, 425-439 (NASA SP-50, Washington, 

D.C., 1964). 

[10] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrody- 
namics of Continuous Media (Pergamon Press, Oxford, 1984) 
[11] For y < 0, B x (x, y, t) is always equal to -B x (x, \y\, t) 



